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A METHOD OF DERIVING FREQUENCY -RESPONSE DATA FOR MOTION 
OF THE CENTER OF GRAVITY FROM DATA MEASURED ON 
AN AIRCRAFT AT LOCATIONS OTHER THAN 
THE CENTER OF GRAVITY 


By John M. Eggleston 


SUMMARY 


A method is presented for deriving time-response and frequency- 
response data for angle of attack and normal acceleration at the center 
of gravity of an aircraft when these data are measured at locations on 
the aircraft other than the center of gravity and when the pitching 
velocity is not measured. The method involves the calculation of transfer 
functions for the measured quantities and operation on these transfer 
functions to derive the transfer function of the desired quantities at 
the center of gravity. Time-response data at the center of gravity are 
then obtained by applying a control input to the corrected transfer func- 
tion. Basic aerodynamic relationships and the forms usually assumed for 
the transfer functions are used in the derivation. A numerical example 
is included to illustrate the application of the method. 

The method appears to have particular application to missiles and 
models of aircraft where the number of telemetering channels available 
may limit the number of quantities that can be recorded or where, because 
of space limitations, instruments cannot be placed in the most desirable 
locations . 


INTRODUCTION 


In the determination of the dynamic stability of aircraft, it is 
often desirable to flight test the aircraft or scale models of the air- 
craft in order to determine transfer functions or frequency responses. 

In the study of longitudinal stability, for example, the transfer func- 
tions relating angle of attack, pitching velocity, and normal acceleration 
to the control deflection may be desired. Frequently, however, facilities 
for recording the data necessary for a general stability analysis are 
limited. Guided missiles and rocket -propelled scale models that require 
telemetering facilities often present such a problem. Furthermore, 
distribution of such items as fuel and propulsion units in either models 
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or full-scale aircraft may preclude the possibility of locating recording 
instruments in the most desirable positions. For these reasons, test con- 
ditions under which a limited number of pertinent quantities are measured 
on the aircraft at positions other than those suitable for stability 
studies are frequently experienced. 

When aircraft motion is confined to two degrees of freedom, any two 
measurements (such as two accelerometers; one in the nose and one in the 
tail) will suffice. It is customary, however, to define the short-period 
longitudinal motion of an aircraft with angle of attack and normal accelera- 
tion measured at the center of gravity. When any two independent quantities 
are measured, time histories of the quantities desired for defining the 
motion may theoretically be calculated from the measured transient response 
data. The procedure, in general, requires differentiating or integrating 
one or more of the recorded quantities and applying corrections point by 
point to the measured curves. This process is laborious, however, and may 
introduce inaccuracies not present in the original data. As an alternative, 
an approach which appears shorter and more accurate would be to calculate 
transfer functions for the measured quantities and to operate on these 
transfer functions to derive the transfer functions and transient time 
responses for the desired quantities. This approach is used in the method 
described in the present paper. 

The derivation of the relationships between angle of attack and 
normal acceleration at the center of gravity and these quantities measured 
ahead, of (or rearward of) the center of gravity is presented for the case 
in which pitching velocity is not known. (The method is greatly simpli- 
fied when pitching velocity is known and consideration of this condition 
is included.) Pitching velocity is calculated from angle -of -attack and 
normal-acceleration transfer functions related to the center of gravity. 

An application of this method is illustrated by use of data obtained 
from drop tests of a freely falling body. 


SYMBOLS 


c 

C L 


c L - 


'Da 


mean aerodynamic chord 
lift coefficient 

dC T 

variation of lift coefficient with angle of attack, 

ox 

variation of lift coefficient with rate of change of angle of 
dCL 

attack, , 

dDo. 


sec 
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Cl D9 

% 

0m 

Cm Da 


variation of lift coefficient with 
variation of lift coefficient with 


pitching velocity, 
elevator deflection. 



d c L 


sec 


pitching-m o ment coefficient 

variation of pitching-moment coefficient with angle of attack, 

Sa 


variation of pitching-moment coefficient with rate of change of 

Sc 

angle of attack, — S, sec 

SDa 






variation of pitching-moment coefficient with pitching velocity, 
3c m 


dDO 


sec 


variation of pitching moment with elevator deflection, -^21 

c)S 


differential operator. 


dt 


g 

Ky 


acceleration due to gravity, ft/sec 2 


nondimens ional radius of gyration of aircraft about axis of 
pitch 


*1 


m 

n 

t 

V 


distance between center of gravity and angle-of -attack vane 
(positive forward), ft 

distance between center of gravity and normal accelerometer 
(positive forward), ft 

numerical integer 

numerical integer or normal acceleration, g units 
time, sec 

true airspeed, ft/sec 


a 


angle of attack, radians 



k 
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6 

9 


flight-path angle, radians 
elevator deflection, radians 
angle of pitch, radians 

relative-density factor, Mass 

Air density x Wing area x c 


phase angle, deg 
angular frequency, radians/sec 
performance operator notation 
Subscripts : 

eg center of gravity 

p location of normal-acceleration pickup 

v location of angle-of -attack vane 

Transfer-function coefficients: 


0 

oo 

[] 


A ° 2(1 fl k Y 2 ( 2 >* f + ° L d 0 ) 

B ’ 2,1 v(v ■ C ®Da) + C Ib9 Cii ‘I}x ■ Cl 0a Cm D9 

0 = f V + 

G - ^ 

H = ° m b V CL De) + C % Cm D3 

J ■ I + C LDa) " < W. 

K - %S> - Vi® 
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N - I S,(%a + c %>) - i + °“De) 

p ■ ifas> - v%) 


E 

F 

X . 

Y 

Z 


coefficients of transfer functions obtained at locations other 
than the center of gravity and corresponding, respectively, 
to coefficients G, H, M, N, and P of transfer functions 
at the center of gravity. 


ANALYSIS 


Consider an aircraft having an angle-of -attack vane located a dis- 
tance Zq ahead of the center of gravity and a normal-acceleration 

sensing device located a distance I 2 ahead of the center of gravity. 

The pitching velocity is not recorded. It is assumed that the aircraft, 
when disturbed from trim, will remain within the linear range of the 
pertinent stability derivatives, that the perturbation angles are small, 
and that the structure is rigid. The sign convention which relates the 
quantities used herein is shown in figure 1. From this figure the 
following basic longitudinal relationship may be Been to exist: 


gn « V D7 = V(D9 - Da) (l) 

The relationship between the normal acceleration at the center of gravity 
and the point at which the normal acceleration is recorded may be written 
as follows: 


-eg - -p - l -i (2> 

A similar relationship may be obtained by relating the angle of attack 
at the center of gravity to the angle of attack measured at the angle-of - 
attack vane 


°*g 



D0 


( 5 ) 
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In cases in which the pitching velocity is known, equations (2) 
and (3) in transfer-function form may he solved to obtain the transfer 
functions of angle of attack and normal acceleration at the center of 
gravity. If the pitching velocity is not known, further transformation 
of the equations is required. Multiplying equation (3) by VD/g gives 

-(Da) = —(Da.) + ^ D 2 6 (k) 

g C S g v g 

and combining equations (l) and (4) to eliminate (Da) C g gives 

"■« + I (Da) - ' (k ' it )* 6 (5) 

The combination of equations (2) and (5) to eliminate D 2 8 gives 


“eg ~ 


1 

(*1 - z 2 )d - V 




( 6 ) 


which is the relationship between the two known quantities of angle of 
attack and normal acceleration located at points other than the center 
of gravity and the unknown quantity of normal acceleration at the center 
of gravity. 

For the relationship of angle of attack at the center of gravity in 
terms of known quantities, equations (2) and (3) may be combined to 
eliminate 9 


a cg “ “v ” 


I 2 VD 


(“eg " ”p) 


(7) 


where and Dp are determined from flight records and r^g is given 

by equation (6) . 

Rearranging the relationships of equations (6) and (j) gives a new 
set of relationships between the known and unknown quantities as 


“eg ' 


rrj 




(8) 


n cg - ”p - (“eg - “v) 


(9) 
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Equations ( 6 ) and (7), or ( 6 ) and ( 8 ), or ( 8 ) and ( 9 ) may be used in the 
determination of angle of attack and normal acceleration at the center of 
gravity inasmuch as these relationships are interchangeable. 


Transfer -Function Concept 


In the application of this method, equations ( 6 ), ( 7 ), ( 8 ), and ( 9 ) 
are used in transfer-function form. The reader is assumed to be familiar 
with the relationship between the time and frequency re alms and with the 
transfer -function concept of expressing the dynamic characteristics of a 
system. Reference 1 presents the theory of this concept and references 2 
and 3 describe several methods of obtaining transfer functions from flight 
data. Equations ( 6 ) and ( 7 ) written in terms of the transfer functions 
relating normal acceleration and angle of attack to elevator deflection 
become, respectively. 




(h D - v )|f] 



( 10 ) 


fal 

a" 

*1« / 

H 

k-*cg 

6. 

lv * 2 VD\ 

5 



( 11 ) 


where brackets denote a performance operator. 

It should be pointed out that no assumption need be made as to the 
form of the transfer functions. Frequency response at the center of 
gravity of the airplane may be obtained from the frequency response meas- 
ured at locations other than the center of gravity by the use of equa- 
tions (10) and (ll) (or equations (8) and ( 9 ) written in operator nota- 
tion) when the substitution of jcs for D is made in the equations. 
However, this approach will not be expanded in this paper. 

In order to illustrate the transf er-function approach, equations (10 ) 
and (ll) are solved by using transfer functions as obtained by stability 
theory when the motions of the aircraft are confined to two degrees of 
freedom. This assumption is made when only the short-period longitudinal 
mode is of interest. If It is desired to assume transfer-function forms 
of higher order than are assumed herein, the method should still be 
applicable . 

The equations of motion for two degrees of freedom (constant air- 
speed and no structural modes of oscillation) are expressed as follows: 
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2 „ s§ %%% - C^DS - C^a - C^Da = 0^8 (12) 

2n |(D0 - Da) - Cl^DS - C^a - C^Da = C ^6 (13) 

Through use of determinants, each of the two dependent variables a 
and D0 may be related to the independent variable 5 in the following 
form: 


a 

5 


GD + H 
AD 2 + ED + C 


(lk) 


D0 = JD + K 
b AD 2 4- ED + C 


( 15 ) 


and by the vise of equation (l), the relationship between n and 6 may 
be shown to be of the form 


n 


5 = 


MD 2 + HD + P 
AD 2 + BD + C 


(16) 


where A, B, C, G, H, and the other transfer-function coefficients 
are composed of comb in ations of the parameters of equations (l P) and ( 13 ) 
and axe defined in terms of these parameters in the list of symbols. It 
may be shown that the denominator of the longitudinal transfer functions 
(the characteristic equation) expresses the damping and frequency of the 
free short-period longitudinal motion of the aircraft and is common to 
all such transfer functions regardless of the location on the aircraft 
at which the variable under consideration is measured or the input and 
response that are related. 

In the application of this method, it is assumed that tr ans fer 
functions relating elevator deflection to the response of angle of attack 
and normal acceleration at points other than the center of gravity have 
been obtained by one of the several means ava ilab le. The transfer func- 
tions discussed herein axe symbolized by the following equations : 
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ED + F 
AD 2 + BD + C 



GD + H 


2 

AD + ED + C 



XD 2 + YD + Z 
AD 2 4 - ED + C 

MD 2 + ED + P 
AD 2 + BD + C 


> 



JD 2 + KD 
AD 2 + BD + C 

J 


(17) 


Reference to equations (2) and (3) shows that the transfer functions at 
the points of measurement and the corresponding transfer functions at the 
center of gravity should he of the same form. 


Transfer Functions at the Center of Gravity 


When the appropriate forms of the transfer functions are used, equa- 
tions (10) and (ll) are expressed in terms of the known and unknown trans- 
fer coefficients. Equation (10), written in terms of equation (17), 
becomes 


MD 2 + ND + P 

1 

- v ) 

XD 2 + YD + Z 

AD 2 + BD + C 

- l 2 )D - Y 

AD 2 + BD + Cj 


Z2VD 2 

ED l + F 

g 

AD 2 + BD + C 


(18) 


Equation (l8) indicates that, with the approach used herein, a 
quadratic numerator and denominator on the left-hand side of the equation 
have been set equal to a cubic numerator and denominator on the right- 
hand side of the equation. Obviously, there must exist a linear factor 
common to both the numerator and denominator of the right-hand side of 
equation (l8) which can be canceled to leave an expression with only 
quadratic factors. 



10 


NACA TN 5021 


Multiplying both sides of equation (l8) by the factor - Z 2 )d " V 

and equating coefficients of equal powers of the operator D yields four 
relationships 


- VP = - VZ (19) 

- Z 2 )P - VN = - VY + IjZ (20) 

^2 2. - 2 2 )n - VM = - VX + IjY + F (21) 

(2 X - 2 2 )m = 2 X X + e (22) 


A second set of relationships may be obtained by following a similar 
process with equation (ll) . Substitution of the transfer-function forms 
of equation (17) into equation (ll) gives 


& 


+ H 

T ED + F 1 Z lS 

(md 2 + 

HD + p) - (xD 2 + YD + z) 

BD + C_ 

[AD 2 + BD + cj 2 2 VD 


AD 2 + BD + C 


(25) 


and equating like coefficients gives the following three relationships : 


*lg 

T£7 


H = F - 


G = E - 


2oV 


(P - Z) 

(24) 

|(N - Y) 

(25) 

!(m - x) 

(26) 


It is immediately obvious that equations (19) and (24) are equivalent 
to a single equation P = Z. Furthermore, examination of equation (22) 
shows that the right-hand side of equation (2 6) may be replaced by the 

or 

expression - S M; therefore 


G = - £ M 
y 


( 27 ) 
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This relationship can also be seen by referring to tbe expressions for 
G and M in tbe section entitled "Symbols." Hence, equations (22) 
and (2 6) are parametric representations of equation (27) and may be 
replaced by this single equation. For purposes of calculation, tbe 
equations which should be used are equations (19)* (20), (21), (25)* 
and (27). Although these equations form the basic linearly independent 
set of equations, equations (22) and (2 6) would also be satisfied if the 


transfer functions 



and 



were known exactly. 


However, in 


the practical application where flight data are used, inaccuracies are 
inherent in the determination of the measured coefficients. In the 
procedure used to determine coefficients relating to motion of the center 
of gravity, the inherent inaccuracies should be minimized. 


Both equation (22) and equation (2 6) contain the coefficient E 
which is the coefficient of the first power of the operator D in the 
angle-of -attack transfer function measured at the vane. Experience has 
shown that E and the associated coefficient of the transfer function at 
the center of gravity G are usually so small that they have little 
effect on airplane response and are therefore usually obtained with low 
percentage accuracy. Hence, although equations (22) and (26) should be 
numerically compatible with the other equations (equations (19) to (27)), 
this numerical check may be poor because of the inherent inaccuracies in 
the measurement of the coefficient E. 


The four relationships, equations (19)* (20), (2l), and (25), may be 
rearranged slightly to express each unknown coefficient as a function of 
the known coefficients. These expressions are 



F + z 

(28) 

V 2 

P = z 

(29) 

Y - — Z 

bo) 

V 

T y - S( l i - *s ) 8 

(31) 


The five equations, equations (27) to (31)* express the coefficients of 
the numerators of the angle-of -attack and normal-acceleration transfer 
functions related to the center of gravity in terns of measured transfer- 
function coefficients, the distance between their point of measurement 
and the center of gravity, and the airspeed at which they were measured. 
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The inverse Laplace transformation may be applied to the determined trans- 
fer functions in order to obtain time-response data at the center of 
gravity in response to control inputs that have Laplace transformations. 


Pitching-Velocity Transfer Function 

The transfer function relating pitching velocity to elevator deflec- 
tion may be obtained in a manner similar to that described for obtaining 
transfer functions involving normal acceleration and angle of attack. 
Starting with any one of the three basic equations involving pitching 
velocity (equations (l), (2), and (3)), the method of equating transfer- 
function coefficients of like powers of the operator D is applied. In 
order to express the coefficients J and K in terms of the measured 
coefficients only, equations (27) to (51) derived in the preceding section 
are then substituted for the derived coefficients to give 

J = F + — Y + ■%(*! - l 2 ) Z (32) 

V V 2 ' 

K = £ Z (33) 

The equations used to determine the unknown coefficients (equations (28) 
to (33)) may be seen to express each unknown coefficient as a function of 
the known coefficient of the same power of D, plus a correction term. 

Other relationships between the coefficients may be readily obtained by 
combining in various ways the equations already presented. The effects 
of using other groups of equations on the accuracy of the numerical results 
are mentioned in the following section. 


EXAMPLE OF APPLICATION 


As an illustration of the method presented herein, use is made of 
data obtained from drop tests of a freely falling model for which the 
response in angle of attack to a step elevator input was measured 5-51 feet 
ahead of the center of gravity (iqj and for which the response in normal 
acceleration was measured 2.165 feet ahead of the center of gravity (^2) • 
The pitching velocity was not recorded. The actual flight data were 
recorded at a constantly decreasing altitude and increasing Mach number. 
These conditions were reflected in the time responses as varying frequency, 
varying damping, and a slight trim change. The response data, however, 
were adjusted for the condition of constant altitude, constant Mach num- 
ber, and zero trim change and were approximated with the following 
analytical expressions for a true velocity of 885 feet per second: 
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<x(t) = 0.01745 - e -1 * l6t (-0.00486 sin 9 - 932 t + 0.01745 cos 9 . 932 t), radians 
n(t) = 0.238 - e “l*l 6 t ( 0 . 3 ^ 6 ! sin 9-932t + 0.17647 cos 9-932t), g units 
6 (t) = - 0 . 00902 , radians 

A plot of these analytical expressions is shown in figure 2. 

Making the Laplace transformation (ref. 1) of these equations and 
dividing the transformations of the first two equations by the transforma- 
tion of the third equation gives transfer functions relating angle of 
attack and normal acceleration to elevator deflection at the points of 
measurement 


"a) 3-109D - 193-40 radian 

- J v Ir + 2 . 32 D + 99.99 radian 



-6.819D 2 + O. 7266 D - 2637-8 g units 
D 2 + 2.32D + 99-99 radian 


These transfer functions are in the normalized form in which all the 
coefficients have been divided by the inertia coefficient A. The sub- 
stitution of D = jo) in these transfer functions produces frequency- 
response relationships that are plotted in figure 3 as amplitude ratio 
and phase angle. The MCA sign convention used in this paper (fig. l) 
defines a positive elevator deflection as one for elevator trailing edge 
down, and therefore, a positive elevator deflection will, in general, 
produce negative static responses. In order to conform with the usual 
practice of plotting frequency-response curves, the phase angles have been 
shifted l80° (zero phase angle at zero frequency). The unusual trend in 
the phase angles shown in figure 3 (a) is caused by the negative coefficient 
of the D 2 term in the numerator of the transfer function relating normal 
acceleration to elevator deflection. 

The application of equations ( 27 ) to ( 31 ) to the coefficients of the 
transfer functions obtained ahead of the center of gravity gives the 
following coefficients corrected to the center of gravity: 


p/a = - 2637.8 
N/A =7-179 
m/a = 6.207 
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g/a = -0.226 
h/a = -194.00 


In order to illustrate the effect on the coefficients M and G 
■when equations involving E are used, calculations of the coefficients 
were made with equations (22) and (26). Equation (22) gave 

' m/A = 44.07 

Equation (26) gave 

G/A = -1.6 (m/a = 44.07) 

g/a = 1.9 (m/a = 6.207) 

These values of the coefficients M and G differ considerably from 
those found by use of equations (27) and (31 ) and give illogical values 
of amplitude ratio at high frequencies in that a much greater value of 
control effectiveness would be required than could be obtained from the 
freely falling model. With reference to the stability parameters that 
define G/A, a positive sign for this coefficient is likewise illogical. 

A check on logical values of the high frequency coefficients M/A and 
G/A was made by noting that from stability theory (see transfer-function 
coefficients in the section on symbols) 


M „ c T £ « 

A 2pgc A 2pc 

By the use of these relationships, values of the parameter Cj^. (based 

on wing area) were found to be in good agreement with theory when equa- 
tions not involving E (i.e., equations (27) and (31)) were used for 
obtaining values of G/A and M/A. 

Additional numerical solutions for the coefficients of the transfer 
functions were obtained by the use of several other sets of equations 
which may be derived from combinations of equations (27) to (33 )• All 
of these sets of equations appeared to be less well conditioned for 
obtaining accurate results than the equations presented because these 
sets of equations frequently expressed the values of unknown coefficients 
in terms of known coefficients of other powers of D to which they are 
not directly related physically. In some cases these other combinations 
of equations yielded numerical answers for the unknown coefficients that 
differed from logical values by as much as a factor of 1000. 
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The coefficients of the numerator of the pitching- velocity transfer 
function may be obtained, by using equations (32) and (33) to give 

J/A = -193. 7^ 

K/A = -95.97 


The complete transfer functions can now be given by the following 
equations : 



-0.226D - 194.00 radian 
D 2 + 2.32D + 99 .99 radian 



= 6.207D 2 + 7.179D - 2637.8 g unit 
D 2 + 2.32D + 99-99 * radian 


D8~ l _ -193 . 7 to - 95-97 radian/ second 

_8J p2 + 2-32D + 99- 99* radian 

The frequency response of the model described by these transfer functions 
is shown in figure 4. 

It may be noted that the significant changes in the transfer functions 
resulting from the transformation to the center of gravity occurred in the 
coefficients of the highest order terms (of the numerators). The negative 
sign of the coefficient G and the positive sign of the coefficient M 
obtained in the transformation axe in agreement with stability theory. 

The phase-angle variation of angle of attack related to elevator deflection 
is such that in the frequency range shown (fig. 4 (b) ), the lag appears to 
approach a constant value of 180°. At higher frequencies, however, where 
the high frequency coefficient becomes dominant, the phase-angle variation 
approaches a constant value of 90° lag. It may also be noted that at the 
natural frequency of the freely falling model, the phase-angle variation 
of the ratio of pitching velocity to elevator deflection is about zero; 
whereas normal-acceleration and angle of attack are lagging elevator deflec- 
tion by about 90° * 

The foregoing results may be further illustrated in the time realm 
by the application of the inverse Laplace transformation to the derived 
transfer functions. For an elevator step input of the same magnitude as 
that used in the teBt, 5(t) = -0.00902 radian, the analytical expres- 
sions as functions of time will be 
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a(t) = 0.0175 - 0.0176e“ 1 * l6t sin(9.932t - 96.6°), degrees 

n(t) = 0.2380 - 0.30592e" 1,l6t sin(9.932t + 84 . 605 °), g units 

D 0 (t) = 0.00866 + 0 . 1752 e _1 * l6t sin( 9 . 932 t - 2 . 83 °), radians/seeond 

These time responses together with the time responses from figure 2 are 
shown in figure 5* 

Inspection of figure 5 indicates that the angle of attack was 
affected very little by the transformation ; whereas normal acceleration 
was changed considerably. The initial positive jump in normal accelera- 
tion at t = 0 experienced by the freely falling model at a point 
2.165 feet ahead of the center of gravity as compared with the negative 
jump at t = 0 experienced at the center of gravity is due to the 
effectively instantaneous change in tail loading as a result of the step 
deflection. The angular acceleration produced by application of the 
tail load resulted in a component of linear acceleration at the point of 
measurement which was larger in magnitude than and opposite in direction 
to the actual linear acceleration produced by the tail load, which is 
indicated by the normal acceleration calculated at the center of gravity. 


CONCLUDING REMARKS 


A method is presented for analyzing longitudinal response data when 
pitching velocity is not known and normal acceleration and angle of 
attack are measured at locations other than the center of gravity. By 
means of this method frequency-response data or transfer functions for 
angle of attack and normal acceleration are obtained. From these trans- 
fer functions the entire time response may be obtained for these quanti- 
ties at the center of gravity. This method also provides a means of 
obtaining pitching velocity when normal acceleration and angle of attack 
are not recorded at the same location in the airplane. Although the 
method as applied used transfer functions obtained from stability theory 
for a system with two degrees of freedom, the analysis was general and 
the same principal should be applicable to systems of higher order or of 
a greater number of degrees of freedom. 

The numerical example presented indicates that care must be used in 
selecting a well conditioned set of equations for the coefficients of the 
transfer function if accurate results are to be obtained. 

This method appears to have particular application to missiles and 
models of aircraft where the number of telemetering channels available 
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may limit the number of quantities that can be recorded or where, because 
of space limitations, instruments cannot be placed in the most desirable 
locations . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 29, 1953* 
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Figure 2.- Time history of elevator input and response in angle of attack 
(measured 5*51 feet ahead of center of gravity) and normal acceleration 
(measured 2.165 feet ahead of center of gravity) of freely falling model. 
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Frequency, go , radians/ sec 

(a) Normal-acceleration response 2.165 feet ahead of center of 
gravity related to elevator deflection. 


Figure 3.- Frequency response of freely falling model as obtained ahead 

of center of gravity. 
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(b) Angle-of-attack response 5. 51 feet ahead of center of gravity- 
related to elevator deflection. 


Figure 3»- Concluded. 


Phase angle, 0, Amplitude ratio, jg 

g units/radian 
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Frequency, co , radians/sec 

(a) Normal-acceleration response at center of gravity related to 

elevator deflection. 

Figure 4.- Frequency response of freely f allin g model at the center 

of gravity. 
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(b) Angle -of-attack response at center of gravity related to 

elevator deflection. 


Figure k.- Continued. 
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Frequency, co, radians/sec 

(c) Pitching velocity related to elevator deflection. 
Figure 4.- Concluded. 
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Figure 5.- Time history of elevator input and response in angle of attack, 
pitching velocity, and normal acceleration of the center of gravity 
as compared with the response obtained at locations ahead of the center 
of gravity of a freely falling model. 
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